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We calculate the equation of state in 2+1 flavor QCD at finite temperature with physical 
strange quark mass and almost physical light quark masses using lattices with temporal extent 
N T = 8. Calculations have been performed with two different improved staggered fermion 
actions, the asqtad and p4 actions. Overall, we find good agreement between results obtained 
with these two 0(a 2 ) improved staggered fermion discretization schemes. A comparison with 
earlier calculations on coarser lattices is performed to quantify systematic errors in current 
studies of the equation of state. We also present results for observables that are sensitive 
to deconfining and chiral aspects of the QCD transition on N T = 6 and 8 lattices. We find 
that deconfinement and chiral symmetry restoration happen in the same narrow temperature 
interval. 

1 Introduction 

A detailed and comprehensive understanding of the thermodynamics of quarks and gluons, e.g. 
of the equation of state is most desirable and of particular importance for the phenomenology of 
relativistic heavy ion collisions. In particular, the interpretation of recent results from RHIC on 
jet quenching, hydrodynamic flow, and charmonium production^ rely on an accurate determi- 
nation of the energy density and pressure as well as an understanding of both the deconfinement 
and chiral transitions. For vanishing chemical p otential , which is appropriate for experiments 
at RHIC and LHC, lattice calculations of the as well as the transition temperature^ 

can be performed with an almost realistic quark mass spectrum. In addition, calculations at 
different values of the lattice cutoff allow for a systematic analysis of discretization errors and 
will soon lead to a controlled continuum extrapolation of the EoS with physical quark masses. 

2 Improved actions and the calculational setup 

Studies of the QCD equation of state are most advanced in lattice regularization schemes that 
use staggered fermions. In this case, improved actions have been developed that reduce 0(a?) 
discretization effects efficiently. Such a reduction is mandatory because a quantitative lattice 
determination of the EoS requires rather coarse lattices as one has to take into account that 
pressure and energy density are dimension 4 operators, which are difficult to calculate on fine 
lattices. At tree level which is relevant for the high temperature phase both the asqtad and the 
p4 actions were found to give rise to only small deviations from the asymptotic ideal gas limit 
already on lattices with temporal extent N T = 6. At N T = 8 the differences from the continuum 
Stefan-Boltzmann value are at the 1 % level At moderate values of the temperature non- 
perturbative effects contribute to the cutoff dependence. In particular in the hadronic phase at 
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Figure 1: The trace anomaly, (e — 3p)/T 4 on N T = 6 and 8 lattices. Results are obtained with the p4 and asqtad 
actions. The right panel shows the low temperature part of the trace anomaly in detail. Here we show fits to the 
data, as well as a comparrison with the resonance gas model (dashed and dash-doted lines). 



low temperature, the breaking of flavor (taste) symmetry, inherent to staggered fermions away 
from the continuum limit, leads to an 0(a 2 ) distortion of the pseudoscalar hadron spectrum 
and may influence the thermodynamics in the confined phase. In order to judge the importance 
of different effects that contribute to the cutoff dependence of thermodynamic observables, we 
have performed calculations with both p4 and asqtad actions which deal with these systematic 
effects in different ways. We extend previous lattice calculations on N T = 6 lattices (asqtad^ 
and to iV T = 8 lattices^. In these studies the bare strange quark mass (m s ) was tuned 

such that the zero temperature kaon mass acquired a constant value of 500 MeV (p4) and a 
slightly larger value for asqtad (570 MeV). The light quark mass was hold at a fixed ratio to 
the strange one of mi = 0.1m s leading to a Goldstone pion mass of about 220 MeV at T = 0. 
The temperature scale was set by using the Sommer parameter r$ = 0.469(7) fm,^ which is 
determined from the shape of the heavy quark potential. 



3 The trace anomaly, pressure and energy density 

The basic thermodynamic quantity most convenient to calculate on the lattice is the trace 
anomaly in units of the fourth power of the temperature O^/T 4 . This is given by the derivative 
of p/T 4 with respect to the temperature, 

8W(r) _e-3p d 4 

T 4 = ~f4~ = T QfWl T > ■ W 

As the pressure is given by the logarithm of the partition function, p/T = V hiZ, the cal- 
culation of the trace anomaly requires only the evaluation of rather simple expectation values. 
In Fig. [H we show results for O^/T 4 obtained with both the asqtad and p4 actions. The new 
-/V T = 8 results^ have been obtained on lattices of size 32 3 x 8 and the additional zero tempera- 
ture calculations, needed to carry out the necessary vacuum subtractions, have been performed 
on 32 4 lattices. The N T = 6 results are taken from Ref. 3 (asqtad) and Ref. 4 (p4), respectively. 

We find that the results with asqtad and p4 formulations are in good agreement. In partic- 
ular, both actions yield consistent results in the low temperature range (see Fig [I] (right)), in 
which 0W/T rises rapidly, and at high temperature, T>300 MeV. This is also the case for the 
cutoff dependence in these two regimes. At intermediate temperatures, 200 MeV<T<300 MeV, 
the two actions show differences and cutoff effects are more pronounced. In the transition re- 
gion the cutoff effects can be well accounted for by a shift of the N T = 6 data towards lower 
temperatures. This reflects the cutoff dependence of the transition temperature and may also 




Figure 2: Energy density and three times the pressure as obtained by p4 and asqtad actions on N T = 8 lattices 
(left) and pressure divided by energy density (p/e) and the square of the velocity of sound (cj) (right), calculated 
on lattices with temporal extent N T = 6, 8. The N T — 6 results for the p4 action are taken from Ref. 3. 



subsume residual cutoff dependencies of the zero temperature observables used to determine the 
temperature scale in the transition region. 

In Fig. [5]we show results for the energy density and three times the pressure. Using Eq. ([I]), 
the pressure is obtained by integrating over the temperature. Crosses with error bars 

indicate the systematic error on the pressure that arises from different integration schemes. The 
energy density (e) is then obtained by combining results for p/T and (e — 3p)/T 4 . Finally, we 
show in in Fig. fright) p/e as function of e. From a spline fit to that data we also calculate the 
square of the velocity of sound, c 2 s . 



4 The QCD transition 

The bulk thermodynamic observables p/T 4 , e/T discussed in the previous sections are sensi- 
tive to the change from hadronic to quark-gluon degrees of freedom that occur during the QCD 
transition; they thus reflect the deconfining features of this transition. In a similar vein, the tem- 
perature dependence of quark number susceptibilities gives information on thermal fluctuations 
of the degrees of freedom that carry a net number of light or strange quarks, i.e., Xq ~ (^g)j 
with N q denoting the net number of quarks carrying the charge q. Quark number susceptibilities 
change rapidly in the transition region as the carriers of charge, strangeness or baryon number 
are heavy hadrons at low temperatures but much lighter quarks at high temperatures. 

Another important aspect of the QCD transition is the spontaneous breaking (restauration) 
of the chiral symmetry. In Fig. fright) we show the quantity A/ s , 

which is defined in terms of light and strange quark chiral condensates evaluated at zero and non- 
zero temperature, respectively. This is an appropriate observable that takes care of the additive 
renormalizations in the chiral condensate by subtracting a fraction, proportional to mi/m s , of 
the strange quark condensate from the light quark condensate. To remove the multiplicative 
renormalization factor we divide this difference at finite temperature by the corresponding zero 
temperature difference, calculated at the same value of the lattice cutoff. 

It is evident from Fig. [3l that A/ S (T) varies rapidly in the same narrow temperature range 
as the bulk thermodynamic observables and in particular the light quark number fluctuations as 
shown in Fig. EJleft). Based on this agreement we conclude that the onset of liberation of light 
quark and gluon degrees of freedom (deconfinement) and chiral symmetry restoration occur in 
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Figure 3: The light quark number susceptibility, calculated on lattices with temporal extent N T = 6 and 8 (left) 
and the subtracted chiral condensate normalized to the corresponding zero temperature value (right). The band 
corresponds to a temperature interval 185 MeV < T < 195 MeV. 



the same temperature range in QCD with almost physical values of the quark masses, i.e., in a 
region of the QCD phase diagram where the transition is not a true phase transition but rather 
a rapid crossover. 

Furthermore, we note that the observed cutoff effects in the chiral condensate can to a large 
extent be absorbed in a common shift of the temperature scale. A global shift of the temperature 
scale used for the N T = 6 data sets by 5 MeV for the p4 and by 7 MeV for the asqtad action 
makes the N T = 6 and 8 data sets coincide almost perfectly. This is similar in magnitude to 
the cutoff dependence observed in (e — 3p)/T and again seems to reflect the cutoff dependence 
of the transition temperature as well as residual cutoff dependencies of the zero temperature 
observables used to determine the temperature scale. 
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